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I. 


Abstract 


The  steady,  zonally  symmetric  motion  of  a  shallow 
incompressible  atmosphere  on  a  rapidly  rotating  earth  subject 
to  an  equator-ward  temperature  gradient  is  studied.  The  assump¬ 
tions  made  allow  the  thermodynamics  to  be  treated  separately 
from  the  motions.  The  turbulence  terms  are  modelled  using  mixing 
length  arguments.  Assuming  the  turbulence  length  scale  is  small 
compared  to  the  earth's  radius  and  the  motion  is  slow  compared  to 
the  earth's  rotation  speed,  the  north-south  geopotential  gradient 
drives  the  eastward  winds  in  the  "interior"  temperate  regions, 
that  is,  not  near  the  surface,  equator  or  poles.  The  meridional 
winds  in  the  interior  are  driven  by  the  turbulence  generated  by 
the  shear  in  the  eastward  winds.  Near  the  equator  the  advective 
terms  become  comparable  with  the  rotational  terms,  but  the 
turbulence  terms  remain  unimportant.  The  motions  there  show 
trade  winds  at  the  equator,  changing  to  eastward  at  a  predicted 
latitude  of  about  15°.  The  meridional  motion  takes  the  form  of 
Hadley  cells  with  rising  at  the  equator  and  sinking  again  at 
about  21°.  The  Hadley  cell  and  the  temperate  region  are 
connected  through  a  vertical  layer  of  turbulence.  Near  the  poles 
the  advective  and  turbulent  terms  become  comparable  with  the 
rotation  terms.  Surface  Ekman  layers  complete  the  picture. 


Preceding  page  blank 
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1.  Introduction 


When  viewed  from  outer  space,  the  earth's  atmosphere,  as 
evidenced  by  the  cloud  patterns,  shows  no  signs  of  symmetry  or 
regularity.  However,  analysis  of  wind  and  temperature  data  over 
long  periods  of  time  shows  that  on  the  average  the  atmosphere 
behaves  in  a  regular  manner.  Prevailing  westerlies,  trade  winds 
and  the  Hadley  circulation  are  gross  features  of  the  motion 
which  are  observed  from  average  data. 

In  1735,  George  Hadley  initiated  the  idea  that  solar 

heating  at  the  equator  forced  the  air  there  to  rise,  and  hence 

he  conjectured  that  the  air,  once  aloft,  travelled  to  the  poles, 

where  it  sank  back  to  the  surface  to  journey  toward  the  equator 

again.  Subsequent  observations  by  Ferrel  around  the  1880's 

suggested  that  the  air  rising  at  the  equator  did  not  travel  to 

the  poles,  but  sank  back  to  the  surface  at  about  30°  north  or 

south  latitude.  In  addition,  the  cold  air  sinking  at  the  poles 

rose  again  at  about  60°  north  or  south  latitude.  In  between 

o  o 

these  cells  was  a  single  cell  rising  at  60  and  sinking  at  30 
in  each  hemisphere.  The  conjectures  mechanism  to  drive  the 
temperature  cell  was  friction  between  the  adjacent  thermally 
driven  cells;  while  the  eastward  winds  therein  were  the  result 
of  the  rotational  (Coriolis)  force  interacting  with  the  meridio¬ 
nal  motions.  Because  of  the  appearance  of  three  cells  in  each 
hemisphere,  mechanisms  to  explain  the  motions  by  explaining  the 
presence  of  each  cell  have  been  called  tricellular  theories. 

The  aforementioned  tricellular  theory  of  driving  the 
motions  has  fallen  into  disfavor,  one  reason  being  that  if  the 
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heating  drives  the  motions  in  the  tropical  and  polar  cells,  the 
winds  in  the  temperate  cell  are  necessarily  weaker  since  it  is 
driven  by  the  others.  But  the  prevailing  westerlies  observed  in 
the  temperate  cell  are  much  too  strong  to  be  secondary  or 
frictionally  forced  winds. 

A  recent  popular  explanation  of  the  driving  mechanism  is 
the  wave  theory  (Pfeffer,  1964).  In  this  theory,  the  north-south 
temperature  gradient  induced  by  solar  heating  drives  wave 
cyclones  and  anticyclones.  The  nonlinear  interactions  of  these 
waves  cause  a  mean  zonal  flow.  The  zonal  flow  in  turn  drives  the 
meridional  cells  through  the  mechanism  of  the  Coriolis  force. 

We  give  a  new  explanation  for  the  general  circulation  in 
which  the  north-south  temperature  gradient  drives  the  zonal 
motions  in  the  temperate  cell  through  the  geopotential  gradient. 
The  zonal  motions  in  the  tropical  and  polar  cells  are  driven  by 
the  zonal  motion  in  the  temperate  cell  through  friction  between 
the  cells.  More  precisely,  the  eastward  motion  at  the  poleward 
edge  of  the  tropical  cell  causes  an  eastward  motion  inside  the 
tropical  cell,  at  least  for  some  distance.  The  Coriolis  effects 
cause  a  strong  meridional  motion  and  the  trade  'in^s.  The 
meridional  motions  in  the  temperate  cell  are  driven  by  frictional 
interactions  with  the  mean  zonal  flow  and  the  tropical  and  polar 
meridional  motions. 

In  order  to  describe  the  motions,  we  shall  use  the  fluid 
mechanical  and  thermodynamical  equations  for  the  zonally  sym¬ 
metric  (i.e.  independent  of  longitude)  flow  of  a  shallow  layer 
of  fluid.  Moreover,  following  Saltzman  (1968),  we  shall  make  e 
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Boussinesq-like  approximation,  neglecting  density  variations 
except  in  the  vertical  momentum  equation.  The  shallow  atmosphere 
assumption  allows  us  to  use  the  hydrostatic  pressure  equation 
for  vertical  momentum  balance.  This  assumption  also  allows  the 
neglect  of  the  vertical  velocity,  except  in  terms  involving 
vertical  derivatives,  which  are  large.  The  hydrostatic  pressure 
equation  allows  us  to  use  the  pressure  as  the  vertical  coordinate, 

A 

replacing  ?,  the  height  above  sea  level.  The  relevant  vertical 
"velocity"  is  oo  =  dp/dt,  the  material  derivative  of  the  pressure. 
We  shall  refer  to  tu  as  the  vertical  p-velocity,  and  to  w  =  dz/dt 
as  the  vertical  z-velocity.  The  pressure  gradient  terms  are 
expressed  in  terms  of  the  geopotential,  defined  as  $  =  gz(J,p), 
where  z((p,p)  is  the  height  of  the  isobar  surface  of  pressure  p. 

The  average  motions  which  we  wish  to  describe  are,  to  a 
large  degree,  both  steady  and  zonally  symmetric.  Thus  we  shall 
assume  that  time  and  longitudinal  derivatives  may  be  neglected. 

When  the  equations  are  appropriately  nondimensionalized, 
there  appear  two  important  parameters,  the  Rossby  number 
Ro  =  U/SOa,  and  the  Ekman  number  £  =  ©/2QU.  Here  U  is  the 
velocity  scale,  O  is  the  rotation  rate  of  the  earth,  a  is  the 
radius  of  the  earth,  and  ©  is  the  scale  of  the  turbulent  flux  of 
momentum  per  unit  mass.  In  the  earth's  atmosphere,  the  Ekman 
number  and  the  Rossby  number  are  both  small,  typical  values  being 

r  _3  _  p 

L  =  10  ,  Ro  ~  10  .  The  presence  of  these  small  parameters  in 

the  problem  allows  the  hope  of  meaningful  approximate  solutions 
valid  in  different  parts  of  the  meridional  ($-p)  plane. 
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For  example,  the  Ekman  number,  which  measures  the  impor¬ 
tance  of  turbulence,  is  qar.te  small,  and  presumably  multiplies 
the  highest  derivatives.  This  implies  that  the  regions  where 
turbulence  is  :.mportant  are  boundary  layers  which  are  thin  com¬ 
pared  to  the  scale  of  global  motions.  The  Rossby  number,  which 
measures  the  importance  of  auvection,  is  also  small,  although 
not  as  small  as  the  Ekman  number.  Thus  the  advection  terms  are 
important  only  in  relatively  small  regions  compared  to  large 
scale  motions.  The  regions  where  advection  is  important,  how¬ 
ever,  are  large  compared  to  the  turbulent  or  Ekman  layers. 


2.  Equations  of  Motion 

The  equations  of  motion  and  heat  balance  for  the  zonally 
symmetric  flow  of  a  shallow  atmosphere  (Saltzman,  1968)  are 


(2.1) 


afl-  (f .  ...afl  i)0  ■  x ,  0 , 

°P  a  cos  f 


(2.2) 


+  (f  +  i22Ll  G)G  + 

ad*  dp  a  a4 


-  Y  =  0  , 


(2.3) 


^  +  T  T  ==  0  , 
dp  p 


(2.4) 


_  +  dv  c 


.  ^  rs 

dp  a  cos  d<f> 


£s.  =  0 

A  A  ^  > 


(2.3) 


f.  dT  ,  ,  dT  RT  £ 

v  — Y  +  - -  Q  =  0 

adt  dp  cpp 
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A  carat  is  used  to  danota  dimensional  variables.  Hara  u,  v  ara 

A 

tne  eastward  and  northward  components  of  velocity,  (f>  is  the 
latitude,  p  is  the  pressure,  cd  =  dp/dt  (t  is  time),  6  is  the 
gravitational  potential,  a  is  the  planet  radius,  T  is  the  temper¬ 
ature,  R  is  the  gas  constant,  and  Cp  is  the  heat  capacity  of  the 
atmosphere.  Also,  f  -  20  sin  J  is  the  Coriolis  parameter,  where 
O  is  the  angular  velocity  of  planetary  rotation.  The  quantities 
X,  y,  Q,  to  be  specified,  represent  turbulent  diffusion  and 
thermal  forcing  of  the  atmospheric  motions. 

It  is  not  our  goal  to  study  the  thermodynamics  of  the  atmo¬ 
sphere.  We  shall  assume  that  we  can  prescribe  the  gravitational 
potential  $(lj>,  p).  This  assumption,  along  with  the  Boussinesq 
approximation,  effectively  separates  the  thermodynamics  from  the 

/v 

mechanics,  since  by  determining  the  motions  in  terms  of  4>,  we  can 
eliminate  the  velocity  components  from  the  thermodynamic  equa¬ 
tions,  yielding  equations  for  $,  T  and  other  appropriate  vari¬ 
ables.  Saltzman  (1968)  discusses  a  similar  approach. 

The  region  in  space  where  these  equations  govern  the  motion 

is 


0  ^  P  ^  Ps(4>)  * 

A 

where  p  ( (f> )  is  the  pressure  at  the  surface, 
s 

The  boundary  conditions  which  are  needed  to  complete  the 
system  fall  into  two  categories,  turbulent  (viscous '  and  symmetry. 
The  turbulent  boundary  conditions  are 


(2.6a,  b) 


/v 


(2.7) 


u  =  v  =  0  on  p  =  ps(ij))  , 


no  stress  on  p  =  0  . 


Because  the  level  p  =  0  actually  corresponds  to  z  =  oo,  an 
unrealizable  situation,  we  snail  not  attempt  to  satisfy  (2.7). 
This  boundary  condition  can  presumably  be  satisfied  by  inclusion 
of  a  turbulent  boundary  layer  at  p  =  0.  A  more  realistic 
approach  to  the  boundary  condition  at  the  "top"  of  the  atmosphere 
might  be  to  apply  a  condition  at  some  p  =  p  <<  p  (  1 )  corre- 

•I  v> 

sponding  to  a  shallow  atmosphere,  say  at  about  50  km.  The 
appropriate  condition  should  reflect  the  presence  of  the  atmo¬ 
sphere  above  50  km. 

The  symmetry  boundary  conditions  are 


(2.3a,b) 

(2.9a,  b ) 

Let  us  nondimensicnalize  the  equations  with 


v=~  =  o  at  <^  =  0  , 

dtf> 

^  A  I  TT 

v  =  u  =  0  at  <f>  =  ±  J 


(2.10) 


r 


u  =  Uu  , 


v  =  Vv  , 


P  =  P0P  ,  i  =  4>c  +  9 


“  =  (P0V/a^o)oo  , 


$  =  c I>o$  , 


(X,Y)  =  0(X,Y)  , 
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where  U  is  a  velocity  scale,  <j>  is  the  latitudinal  scale  of  the 
motion,  p  is  a  pressure  scale,  taken  to  be  the  surface  pressure 
at  the  equator.  The  scale  of  u>,  P0U/a<f>0,  is  suggested  by  the 
continuity  equation  (4).  We  also  assume  that  is  an  appropri¬ 
ate  scale  for  the  gravitational  potential,  and  0  is  a  sca,e  for 
the  turbulent  stresses. 

We  shall  denote  0/20  U  by  £  ,  and  U/20a  by  Ro.  We  call 
these  two  parameters  the  Ekman  number  and  the  Rossby  number, 
respectively.  It  is  important  to  note  that  the  Rossby  number 
defined  h^re  is  not  that  used  by  Charney  (1948,1963).  He  includes 
the  quantity  sin  4>  in  his  Rossby  number.  The  Rossby  number  we 

A 

have  defined  is  a  constant,  independent  of  (j>. 

In  terms  of  these  parameters  the  nondimensionalized  equa¬ 


tions  become 


(2.11) 


du  cos  ( (j>c  +  4>0<1>) 
cos(<t>c  +  4>04> ) 


-  a  sin(  (J>c  +  (}>0<t>)v  =  (fx  , 


(2.12)  (v  ||  +  a>  ||)  +Ro  tan(<|>c  +  4>04> )u2  +  sin(<j>c  +  i>0+)u 

*o  ^ 

=  lfY  '  (4>c  +  M'p)  > 

^  av  cosU,c^on  _  n 

"5p  cos  ( <f>c  +  4>0<|>)d<i> 


3 .  Turbulence  Model 

Let  us  now  consider  models  for  the  turbulence  terms.  The 
northward  transport  of  eastward  momentum  due  to  turbulence  is 
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(3.1) 


1 


/N 


( (v'u* ) 


cos 


a  cos 


|  d(j> 


t) 


) 


where  a  prime  denotes  the  fluctuation  from  the  zonally  averaged 
flow,  and  the  bar  denotes  the  zonal  average.  We  shall  assume  a 
mixing  length  hypothesis  based  on  the  assumption  that  differences 
in  the  eastward  velocity  are  most  effective  in  supporting  the 
turbulence.  For  a  discussion  of  mixing  length  theory,  see  Hinze 
(1959). 

We  assume  that 


(3.2a) 

and 

(3.2b) 


u' 


h  du 
a  d$ 


where  and  are  mixing  lengths.  The  northward  transport  of 
eastward  momentum  due  to  turbulence  is 


(3.3) 


2  2 

LHU  1 


N 


a "’dr'  cos 
o 


_d_  [•  2  du  du 

|  d<t>  1  «  ^  ^ 


COS  (|)} 


where  is  the  horizontal  size  scale  of  the  eddies,  and  AH  is 
an  average  dimensionless  mixing  length,  defined  by 


(3.4)  L^jj  -  J  |  . 

The  appearance  of  the  minus  sign  is  accounted  for  by  the 

following  argument.  If  du/d<|>  >  0  and  >  0,  then  u’  >  0;  but 

the  average  momentum  transport  is  southward,  since  we  are  moving 

Sj  *|  A 

a  parcel  with  excess  momentum  - ^  southward,  so  that  v*  <  0. 

A  A 

Similar  arguments  for  du/d^  <  0  lead  to  (3.3). 
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The  vertical  transport  of  eastward  momentum  is 


(3.5) 


If  we  again  use  mixing  lengths,  we  have 


(3.6a) 

and 


u' 


du 

Sp 


(3.6b) 


~  ^  ^u 

“  =  PS  %  -x 

dp 


y 


where  we  now  assume  that  the  velocity  differences  which  are 


effective  in  the  vertical  momentum  transport  are  vertical  differ¬ 
ences,  so  that  the  relevant  derivative  is  du/dp.  Here  £  and 
are  effective  pressure  differences  which  scale  the  vertical 
turbulent  eddies.  The  quantity  pg  converts  vertical  z-velocity 
to  vertical  p-velocity .  This  hypothesis  leads  to  the  following 
expression  for  the  /rrtical  transport  of  eastward  momentum  due  to 
turbulence : 


(3.7) 


*v  - 


2  2 


b 

^P 


tp*S‘$p>8> 


Here  ly  is  the  vertical  size  scale  of  the  eddies  and  A^  is  an 
average  dimensionless  mixing  length,  pQ  is  the  density  at  the 

A 

surface,  and  p  =  pQp. 

For  the  turbulent  transport  of  northward  momentum,  we  shall 
use  similar  considerations.  Again  we  shall  assume  that  differ¬ 
ences  in  eastward  velocity  are  most  effective  in  supporting  the 


turbulence.  This  leads  to  the  expression 


(3.8) 

and 

(5.9) 


7%  Jj 


*0 


where  [xu  and  a...  are  average  dimensionless  mixing  lengths.  We 

n  V 

shall  assume  that  A^,  Ay,  and  Uy  are  constants.  Then,  without 
loss  of  generality,  we  can  take  A^  =  1. 

The  corr^t  scale  0  for  the  turbulence  terms  is  the  larger 
of  the  vertical  and  horizontal  scales. 


(3.10) 


0  =  max 


2  2  2  2 

WU 


By  talking  =  100  km,  and  Ly  =  100  mb,  corresponding  to  a 
vertical  eddy  scale  of  a  couple  of  kilometers,  we  find  that  the 
horizontal  and  vertical  turbulence  terms  are  comparable  if  (f>rj  =  l. 
If  <t>  is  smaller,  the  horizontal  terms  are  dominant.  Unless 
otherwise  specified,  we  shall  use  the  horizontal  turbulence  terms 
for  scaling  purposes.  We  denote  L^/a  by  e.  We  now  assume  that 

C  2  R 

e  <<  1  and  Ro  «  1.  Since  G  =  ^.  .  _2  ,  the  turbulence  terms  will 

<t>0 

be  negligible  compared  to  the  advective  terms  as  long  as  4)c,  >>  e» 
Thus  for  latitudinal  scales  larger  than  a  typical  large  scale 
eddy,  we  can  neglect  the  turbulence  terms.  More  precisely,  the 


role  of  turbulence  has  been  relegated  to  various  boundary  layers 


in  the  flow. 


—  - - —  ■  ■ 

4.  The  Interior  or  Temperate  Approximation 


In  seeking  approximate  solutions,  it  is  natural  to  assume 

4,  =  i,  d>  =  o,  and  a  =  1  and  neglect  terms  of  order  Ro.  This 

Tq  5  Tc 

yields  the  approximation 


(4.1) 

(4.2) 

(4.3) 


-sin  =  0  , 


sin  <|>u  =  -  d$/d<}> 

+  dv  cos  o 
"5p  cos  (J)  d<{> 


0 


Here  we  have  set  =  20Ua. 

Equation  (4.1)  suggests  not  that  v  =  0,  but  that  the  first 

term  in  the  expansion  of  v  is  small  compared  to  1.  An  Ekman 

layer  exDansion  for  the  surface  boundary  layer  (Pedlosky,  1969) 

.  ,  rl/2 

suggests  that  the  next  term  in  the  expansion  is  of  order  £ 
Pedlosky  ignores  the  inertia  terms  and  uses  a  viscous  term  to 
model  the  effects  of  turbulence.  Inclusion  of  the  inertia  terms 
does  not  change  the  result.  However,  at  this  order  in  the  inner 
equation  v  is  still  zero. 

The  appropriate  balance  in  the  inner  region  is  a  =  C  ,  so 

that 

A 

(4.4)  -sin  <f)v  =  X  . 


From  equations  (4.1 -4.4),  the  mechanism  driving  the  motion 
in  the  interior  (temperate)  region  is  as  follows:  the  geo¬ 
potential  drives  the  eastward  motion  through  the  Coriolis  force. 
The  meridional  motions  v  and  00  are  driven  by  friction;  more 
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specifically,  the  turbulent  stress  in  the  eastward  direction 
drives  the  northward  motion.  The  vertical  motion  is  determined 
through  the  continuity  equation. 

In  order  to  determine  the  motions  further,  it  is  necessary 
to  have  expressions  for  the  geopotential  d>  and  the  density  p, 
which  appears  in  the  eastward  turbulent  stress  force  X.  Thermo¬ 
dynamic  considerations  are  important  in  each  expression. 

Since  the  velocity  scale  U  can  be  taken  to  be  the  largest 
value  assumed  by  the  magnitude  of  the  velocity  vector,  it  is 
approximately  (in  dimensional  variables) 


(4.5) 


U  =  max 


20a  sin  <£ 


If  wa  take  values  of  $  from  Schutz  and  Gates  (1971)  the  velocity 
scale  derived  from  the  model  is  about  30  m/sec,  which  gives  a 
Rossby  number  of  about  0.03.  The  eastward  velocity  in  the 
northern  hemisphere,  as  calculated  from  the  data  of  Schutz  and 
gates,  is  shown  in  Figure  1. 

From  the  meridional  motions,  using  the  mixing  length 
hypothesis  of  Section  3,  we  have 


(4.6) 


v  = 


f  1  d  T  Su|  Sul 
n  h  l^T"^  l^rl^|cos 


sin  <{> 


T  "5$  L§?l5?l 


*]+Bi  <!p>2}  - 


where  p 


p  p  -z. 

Vva  P0g 

PoLH 


is  assumed  to  be  of  order  1, 


5.  The  Temperate  Surface  Boundary  Layer 


There  are  several  boundary  layers  needed  to  complete  the 
picture.  It  is  well  known  that  an  Ekman  (turbulent)  boundary 
layer  is  needed  to  satisfy  the  boundary  conditions  at  the  earth's 
surface  (2.6a, b).  In  the  temperate  region,  where  the  advective 
terms  are  negligible,  the  analysis  of  Pedlosky  is  valid. 

Pedlosky,  however,  uses  Newtonian  viscosity  cerms  to  model  the 
turbulence.  We  feel  that  the  effect  of  turbulence  car.  be  better 
modelled  by  the  mixing  length  considerations  introduced  before. 

In  the  temperate  region  we  wish  to  balance  the  turbulence 
terms  with  the  rotation  terms.  This  leads  to  a  pressure  scale  of 
(pogL^  311(1  the  following  equations  of  motion. 


(5.D 


ip  =  ~sin 


<J)V 


(5.2) 

(5.5) 


2  a 
uv  Ip 


((|H)2)  =  sin  iu  +  ^  (i,p8(i)) 
y  d<{> 


> 


dv  cos  (j) 

- x - *■ 


cos  c|>  d<}> 


0  , 


where  Ay  was  taken  to  be  1. 

Since  the  surface  boundary  layer  does  not  shed  any  light 
on  the  global  circulation,  we  shall  not  discuss  it  further  here. 


6.  The  Equatorial  Cell 

The  approximation  <j>  =  1  is  not  valid  near  the  equator, 

A  A 

where  is  small,  or  near  the  poles,  where  tan  <J>  is  large.  To 
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find  solutions  valid  in  these  regions,  we  must  re-scale  the  lati¬ 
tude  by  choosing  <j>0  to  be  small.  The  result  is  an  advective 
"boundary  layer",  which  is  small  compared  to  interior  sizes,  but 
large  compared  to  the  Ekman  layers. 

To  obtain  an  approximation  valid  near  the  equator,  we  let 
(f>  =  (Ro)1/2,  =  0,  und  a  =  1.  The  equations  then  become 


(6.1) 


(6.2) 


(6.3) 


$  +  H  ■  0(r°1/2)  * 

Ro  (v  |f  +  “  H  +  +  If  *_0)  =  0(Ro3/S)  ' 

Ro  ( v  If  +  m  tp  ‘  ^  =  o(R°5//2)  • 


1/2 

In  this  case  we  scale  the  geopotential  by  <J>  =  (Ro)  '  (2Uf)). 


Thus,  to  lowest  order  in  Ro  we  have 


(6.4) 


dv  ,  dcs  _  n 

35  '  ’ 


(6.5) 


du  ,  du 
v  +  “  =  *v  ' 


(6.6) 


dv  ,  dv  1  C*I> 

v  Sf  +  ’  ‘*U  '  8?  *.0 


In  terms  of  the  equatorial  dimensionless  variables,  the  region 
on  which  we  should  consider  (6.4 -6.6)  is 


(6.7) 


-co  <  4>  <  00  ,  0  _<  p  1 


Let  us  attempt  a  solution  to  the  equations  of  motion  (6.4 
6.6)  under  the  conditions  corresponding  to  "spring"  or  "fall". 
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that  is,  symmetric  forcing, 
and  the  equations  become 


If  <J>  is  symmetric,  then  d$/d<f) 


<j>=0 


=  0, 


(6.8) 

v  If  +  “  If  =  - 

(6.9) 

v  1? + «  If  -  • 

(6.io) 

• 

o 

ii 

+ 

>he~ 

The  relevant 

boundary  conditions  are  (2. 

surface  boundary  1  i.yer  is  needed  to  satisfy  (2.6a, b). 

The  balance  in  the  equatorial  region,  as  expressed  by  equa¬ 
tions  (6.8-6.10)  is  the  balance  between  rotation  and  advection. 


The  driving  mechanism  (source  of  energy)  is  the  boundary  layer 
surrounding  the  equatorial  region,  with  the  important  driving 
being  through  a  turbulent  layer  between  the  equatorial  region 
and  the  interior  or  temperate  region. 

The  continuity  equation  (6.10)  implies  that  there  exists  a 
stream  function  ip  with 


Since  we  assume  that  there  is  no  flow  through  the  top  (p  =  0)  or 
the  bottom  (p  =  1)  of  the  atmosphere,  the  lines  p  =  0  or  1  are 
streamlines. 

The  system  of  partial  differential  equations  (6.8-6.10)  is 
hyperbolic;  hence  there  exist  characteristic  curves.  The  equa¬ 
tions  for  the  characteristics  are 
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(6.12) 

d<f> 

dt 

(6.13) 

dp  _ 
dt 

(6.14) 

du 

eft  “ 

(6.15) 

dv 

where  t  is  a  variable  along  a  characteristic. 

We  observe  from  this  system  that  dtp/dt  =  0  so  thac  the 

characteristics  for  the  system  are  the  streamlines. 

From  the  two  momentum  equations  we  see  that  the  horizontal 

2,2 

kinetic  energy  E  =  - — 7^—  is  constant  on  each  characteristic. 

We  shall  use  E  to  parameterize  the  characteristics.  This  implies 
that  Tp  =  ip  ( E ) . 

The  eastward  momentum  equation  (6.14)  along  with  (6.12) 
give  u  =  i^-  +  uo(E).  We  then  obtain  v  from  u2  +  v2  =  2E,  giving 


(6.16) 


2E-^-uoi2-U2 


The  solutions  to  the  <{>  equation  (6.16)  depend  on  the  sign 
of  uQ.  In  order  to  get  north- south  symmetry,  we  must  choose 
u  <0.  This  gives  v  =  0  at  =  0,  and  hence  allows  a  conver¬ 
gence  zone  at  the  equator.  The  resulting  solution  for  the  <{> 
equation  (6.16)  is 

$  t 

(6. .’7)  <Mt)  »  ±<i>0vE)  dn  (~g-,  k) 

where  dn  is  the  delta  amplitude  function,  and  <J>Q  is  related  to 

uo  by 


\ 
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ir  +  uo  =2E' 


(6.18  1 


It  remains  to  determine  the  location  of  the  streamlines. 


To  do  tnis  we  write 


(6.19) 


for  each  <j)  and 


t' 

-I 


dp  +pr 


p0<*> 


i  -  -  la^asi  (8E/aP) . 


We  shall  assume  that  at  p  =  p  (<{)),  v  =  0.  Thus 


(6.20) 


where 


(6.21) 


dtf'/dE' 


E°(+)  ±isE' 

2E0(4.)  -  (u0(E0(«,))  . 


<JE’  +  P„(4>) 


We  now  impose  the  conditions  that  p  =  0  and  p  =  1  are 
characteristics.  We  note  that  since  (f)  =  0  lies  on  both  these 
curves,  the  value  of  E  corresponding  to  these  curves  is  given  by 


(6.22) 


This  yields 


2EX  =  (uo(E1))2  . 


(6.2^) 


P0(<M  =  1/2 


W4 


(6.24) 


This  integral  equation  gives  a  relation  between  ip(E)  and  uq(E). 
Thus  for  this  system,  on  the  domain  given  previously,  there  is 
one  arbitrary  function,  which  we  may  take  to  be  u  (E). 

The  appearance  of  an  undetermined  function  in  the  solution 
is  a  common  occurrence  in  inviscid  flows.  Presumably  the  func¬ 
tion  U0(E)  can  oe  determined  by  including  dissipative  effects, 
as  in  Batchelor  (1956).  The  condition  which  uq(E)  must  satisfy 


Xdt  =  0  , 


E=const 


where  X  -  •“  •  This  leads  to  the  following  functional 


differential  equation  for  uq: 


(6.25) 


Auo  +  Buouo  +  CuoS  +  +  F  =  0  , 


where  a  prime  denotes  differentiation  with  respect  to  E,  and 

A  i<jo2dt/ty?2  , 

B  =^(j)  a?  dt  , 

C  =^axiyty'2  +^aj5dt^"/V/f  -^)fco2dt , 

D  =  -  ^  ^dt/^'  , 


are  functionals  of  u^.  Since  the  equation  is  second  order,  we 
need  two  conditions  on  u  (E1).  These  will  be  discus' ed  in 
Section  7* 


We  shall  not  attempt  to  solve  equation  (6.25)  here. 
Instead,  let  us  examine  the  motions  obtained  by  taking  u  con¬ 
stant.  One  implication  of  this  assumption  is  that  the  eastward 
velocity  is  then  independent  of  the  vertical  coordinate.  This 
assumption  is  somewhat  unreal  physically,  but  allows  us  to  solve 


the  integral  equation  (6.24). 

u 

With  uq  taken  to  be  a  constant,  E1  =  and 

<«0+£a 

Eq(4>)  =  - g - .  The  equation  itself  is  an  Abel's  equation, 

and  its  solution  is  easily  determined  to  be 


(6-26>  &-$-=+=• 

K-se 

Substituting  this  result  in  equation  (6.20)  for  p  gives 


The  resulting  motions  are  shown  in  Figures  2  and  5. 

The  meridional  motions  show  Hadley  cells,  as  observed  in 
the  tropical  regions  of  the  earth's  atmosphere.  These  cells 
extend  to  latitude  <j>  =  ±2(Ro|uo|)^^  in  either  hemisphere.  If 
we  assume  that  the  appropriate  value  of  uq  is  tha„  which  gives 
the  global  maximum  value  of  u  at  the  edge  of  the  cell,  we  have 
uo  "  and  we  find  that  the  cells  extend  to  about  21°  north  or 
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south  latitude.  With  this  assumption,  the  change  from  westerly 
winds  to  trade  winds  occurs  at  <$>  =  ±  /2  (Ro|u  |)  '  ,  corre¬ 
sponding  to  about  15°.  The  predicted  values  agree  reasonably 
well  with  observations. 

Using  considerations  similar  to  those  in  Section  5,  we  find 
that  the  surface  boundary  layer  in  the  tropical  Hadley  cell  is 
governed  by  the  equations 


(6*P8)  jp  =  4  V  +  03  ' 

(6.29)  <fm  +  v  ||  +  a)  ||  , 

|  +  t  -  °  * 


The  vertical  length  scale  for  this  layer  is  (pQgL^-^  (Ro)1/2)1^. 
We  shall  not  discuss  this  layer  further  here. 


7.  The  Turbulent  Connecting  Layer 

From  equations  (6.17)  and  (6.18)  we  see  that  the  solutions 
in  the  Hadley  cell  region  are  defined  only  for 

A  ]/2 

|(J>|  j<  2(Ro|uo(E1 )  |  )  ,  even  when  uq  is  not  assumed  to  be  a  con¬ 

stant.  It  is  this  failure  of  the  solution  which  indicates  that 
a  turbulent  boundary  layer  is  needed  to  connect  the  Hadley  cell 
to  the  temperate  region.  Moreover,  since  equations  (6.8-6.10) 
have  no  forcing  term,  the  driving  mechanism  for  the  Hadley  cell 
must  be  movement  of  part  of  its  boundary.  We  claim  that  the 


Hadley  circulation  is  driven  by  the  turbulent  layer  at  the  pole- 
ward  edges  of  the  cells. 

The  correct  scaling  for  the  turbulent  layer  between  the 
r  >p  cal  Hadley  cell  and  the  temperate  region  is  =  2(Ro)’^ 

and  4)q  =  c.  We  assume  that  e  <  (Ro)1/2.  The  approximate  equa¬ 
tions  become 


(7.1) 

(7.2) 

(7.3) 


& 


f#2)  - 


v 


du 


+  0) 


du 
3p  ' 


U 


2  d  /du,2 
H  "§(£  "cjif 


bv 
3p  ' 


dv  (  dcu 

d4>  "5p 


0 


to  order  (Ro)1/2. 

These  equations  are  difficult  to  solve.  We  shall  satisfy 
ourselves  by  assuming  that  u  is  continuous,  and  noting  that  if 
we  integrate  (7.1)  or  (7.2)  from  0  to  1  in  p,  and  from  -co  to  co 
in  we  derive  the  result  that  across  the  boundary  layer,  the 
quantity 
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is  conserved.  Physically,  the  second  condition  states  that  the 
total  northward  turbulent  flux  of  eastward  momentum  through  the 
layer  is  a  constant. 

Wc  can  estimate  the  thickness  of  the  turbulent  layer  by 
noting  that  its  latitudinal  thickness  scale  is  e,  the  ratio  of 
the  eddy  size  to  the  earth's  radius.  If 


we  assume  an  eddy  size 


of  a  few  hundred  kilometers,  then  e  is  comparable  to  Ro,  the 
Rossby  number,  and  <f) .  is  about  2°.  Thus  we  might  expect  its 
influence  to  be  felt  for  about  10  or  so* 


8.  The  Polar  Region 


The  mechanics  in  the  polar  region  is  somewhat  more  com¬ 
plicated  than  that  of  the  tropical  regions.  If  we  attempt  a 
solution  analogous  to  the  tropical  aavective  boundary  layer 
surrounded  by  turbulent  layers,  the  turbulent  layers  are  as 
large  as  the  advective  region.  This  suggests  that  the  equations 
appropriate  in  the  polar  region  entail  a  balance  between  rota¬ 
tion,  advection  and  turbulence.  In  order  to  model  this,  we 
assume  that  <f>c  -  |,  <|>0  =  Ro,  a  =  1  and  e  =  Ro.  The  resulting 
approximate  equations  are 


(8.1) 

(8.2) 


(8.3) 


dv  .  du>  .  v 

af  +  3p  +  T 


to  order  Ro. 

In  terms  of  the  polar  dimensionless  variables,  the  region 

under  consideration  is  0  >  <j>,  0  <  p  <  Pn  where  pn  -  Ps(^)/Ps(°)« 
In  order  that  the  geopotential  be  smooth  at  the  pole,  we 

must  require  ^t>(^>  P )  /d<t>  = 
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The  equations  valid  in  the  polar  region  are  not  readily 
treatable  by  classical  means.  We  can,  however,  estimate  the 
extent  of  this  polar  region,  which  we  assume  to  be  a  cell,  as  in 
the  tropical  case.  Since  <f>  =Ro,  corresponding  to  a  few  degrees 
latitude,  the  polar  region  should  be  several,  perhaps  10  degrees 
wide.  This  is  somewhat  smaller  than  observations  indicate  for 
the  earth. 


9.  The  Energy  Budget 

A  description  of  the  motions  of  the  atmosphere  would  not  be 
complete  without  an  explanation  of  the  energy  budget.  The  asymp¬ 
totic  approach  taken  here  suggests  that  the  following  process  is 
responsible  for  the  dissipation  of  the  potential  energy  caused 
by  solar  heating.  That  potential  energy  is  in  the  form  of  a 
north-south  geopotential  gradient.  This  geopotential  gradient 
drives  the  eastward  flow  in  the  middle  latitudes,  and  hence  is 
changed  to  the  kinetic  energy  of  the  mean  zonal  flow  there. 
Through  turbulence  generated  by  the  zonal  flow  in  the  form  of 
large  scale  eddies  and  their  related  vertical  overturnings,  the 
kinetic  energy  of  the  zonal  flow  is  changed  into  kinetic  energy 
of  the  meridional  flow. 

The  zonal  flow  in  the  middle  latitudes  also  drives  large 
scale  eddies  at  the  interface  between  the  tropical  Hadley  cell 
and  the  temperate  region.  The  large  scale  eddies  in  turn  drive 
the  motions  in  the  Hadley  cell,  where  the  zonal  motion  is  coupled 
to  the  meridional  motion  through  the  Coriolis  forces. 
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The  zonal  flow  also  drives  the  motions  in  the  polar  cells, 
which  in  general  are  weaker  than  those  in  the  tropical  cell.  For 
a  schematic  representation  of  the  energy  cycle,  see  Figure  4 . 


10.  Jupiter's  Atmosphere 

It  is  of  interest  to  apply  some  of  the  ideas  presented  in 
this  paper  to  the  motions  of  the  atmosphere  of  Jupiter. 

Jupiter's  atmosphere,  with  the  major  exception  of  the  Great  Red 

Spot,  appears  to  be  zonally  symmetric.  In  addition,  the  estima- 

-8  -2 

ted  Ekman  and  Rossby  numbers  for  Jupiter  are  10  and  10 
respectively.  In  that  case,  the  analysis  presented  here  should 
be  able  to  predict  some  features  of  the  motion,  provided  that  the 
atmosphere  meets  the  other  requirements  herein,  such  as  the 
shallow  atmosphere  assumption. 

Since  a  knowledge  of  the  geopotential  is  required  for  pre¬ 
diction  of  the  winds  in  the  temperate  region,  we  shall  not  concern 
ourselves  with  that  part  of  the  atmosphere.  Some  features  of  the 
motions  in  the  tropical  zones,  on  the  other  hand,  can  be  discerned 
from  the  value  cf  the  Rossby  number.  In  particular,  the  motions 
show  a  Hadley  cell,  of  width  approximately  10°.  This  width 
agrees  extremely  well  with  the  width  of  the  equatorial  zone  on 
Jupiter.  The  present  model  also  predicts  trade  winds,  changing 
to  prevailing  westerlies  at  about  8°.  This  prediction  is  not 
observed;  in  fact  the  east-west  circulation  is  observed  to  be 
eastward. 
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11.  Summary 


Let  us  summarize  the  asymptotic  scheme  we  have  presented 
for  appi oximat ing  the  motions  of  the  earth's  atmosphere.  We 
separated  thermodynamic  considerations  from  the  equations  of 
motion  by  assuming  that  the  temperature  field  is  known.  This 
allowed  us  to  treat  the  fluid  dynamical  equations  as  forced  by 
thf  geopotential. 

We  assumed  a  mixing  length  theory  for  the  turbulence  terms. 
By  assuming  that  the  eddy  size  is  small  compared  to  the  earth's 
radius,  we  found  that  the  effect  of  turbulence  is  confined  to 

relatively  small  boundary  layers. 

We  also  assumed  that  the  motions  we  seek  are  slow  compared 
to  the  speed  of  the  earth's  rotation  speed.  This  implies  that 
the  inertial  terms  are  negligible  except  in  layers  which  are  rela¬ 
tively  narrow  compared  to  the  whole  atmosphere,  but  large  compared 
to  the  turbulent  layers.  These  advective  boundary  layers  occur 
near  the  equator  and  the  poles,  and  take  the  form  of  meridional 
cells. 

The  dynamics  of  the  motion  can  be  summarized  as  follows. 

In  the  temperate  region,  the  north-south  geopotential  gradient 
drives  the  prevailing  westerlies.  Turbulent  eddies  drive  the 
meridional  motions.  The  eastward  flow  at  the  temperate  edges  of 
the  polar  and  tropical  cellr  drives  the  motions  in  these  cells. 

In  the  equatorial  cell,  the  balance  between  the  Coriolis  force 
and  the  advection  causes  the  prevailing  westerlies  to  change  to 
trade  winds  near  the  equator,  and  also  creates  the  meridional 
cell  pattern. 
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There  is  a  turbulent  layer  connecting  the  equatorial  cell 
to  the  temperate  zone. 

The  mechanics  in  the  polar  region  is  complicated  by  the 
fact  that  the  turbulent  terms,  the  Coriolis  terms  and  the  advec- 
tive  terms  are  all  of  equal  importance. 

Approximate  solutions  for  the  motion  are  not  found  in 
regions  where  turbulence  is  important.  Presumably  numerical 
techniques  are  needed  for  these  regions. 

To  have  a  complete  description  of  the  atmospheric  dynamics, 
we  need  to  be  able  to  derive  expressions  for  geopotential  £  and 
the  density  p,  which  can  be  expressed  in  terms  of  the  temperature 
x  as  p  =  p/RT.  The  geopotential  and  temperature  are  determined 
througn  the  energy  equation  and  the  vertical  momentum  equation, 
along  with  other  balance  laws  for  quantities  which  are  important 
in  thermodynamic  considerations,  such  as  water  vapor,  cloud 
cover,  etc.  In  general,  these  equations  involve  the  velocity 
components,  for  which  we  have  approximate  expressions  involving 
the  geopotential  and  the  temperature.  In  theory,  we  need  only 
substitute  the  appropriate  expressions  for  the  velocity  components 
where  they  appear  to  yield  a  thermodynamic  system  involving  only 
thermodynamic  variables,  independent  of  the  velocity  components. 

In  fact,  however,  the  expressions  for  the  velocity  components  are 
not  simple  expressions  involving  the  thermodynamic  variables,  and 
so  this  procedure  does  not  lead  to  easily  tractible  thermodyna lie 
equations.  We  shall  leave  the  thermodynamic  considerations  foi 
the  future. 
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An  approximate  description  of  atmospheric  dynamics  could 
be  useful  in  predicting  climate  changes  due  to  various  effects, 
for  example,  changes  in  cloud  cover,  solar  absorption,  et^. 

Most  of  these  changes  alter  the  thermodynamics  of  the  atmosphere. 
In  addition,  the  thermodynamics  is  affected  by  the  motion  of  the 
atmosphere.  What  we  have  presented  is  an  approximate  description 
of  the  motions  derived  from  given  thermodynamic  conditions.  This 
reduces  the  problem  to  consideration  of  thermodynamic  equations. 
Thus  it  is  felt  that  this  work  represents  a  step  toward  a  useful 
analytical  approach  to  global  climatology. 
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